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The generator coordinate method (GCM) casts the wavefunction as an integral over a weighted 
set of non-orthogonal single determinantal states. In principle this representation can be used like 
the configuration interaction (Cl) or shell model to systematically improve the approximate wave- 
function towards an exact solution. In practice applications have generally been limited to systems 
with less than three degrees of freedom. This bottleneck is directly linked to the exponential compu¬ 
tational expense associated with the numerical projection of broken symmetry Hartree-Fock (HF) 
or Hartree-Fock-Bogoliubov (HFB) wavefunctions and to the use of a variational rather than a 
bi-variational expression for the energy. We circumvent these issues by choosing a hole-particle 
representation for the generator and applying algebraic symmetry projection, via the use of ten¬ 
sor operators and the invariant mean (operator average). The resulting GCM formulation can be 
mapped directly to the coupled cluster (CC) approach, leading to a significantly more efficient 
approach than the conventional GCM route. 


I. THE GENERATOR COORDINATE 
METHOD (GCM) 

The generator coordinate method (GCM), initially de¬ 
veloped to treat large amplitude collective motion in nu¬ 
clei pHl] (see also Refs. US) , represents the wavefunc¬ 
tion as an integral over a weighted set of states 

|*GCM> = ydZ/(Z)|*(Z)>, (1) 

each state |<J>(Z)) parameterized by a continuous variable 
Z, the so-called generator coordinate. For bosonic modes 
this structural form is related to the coherent state rep¬ 
resentation Q. Unlike the configuration interaction (Cl) 
or shell model 

I’l'ci) = X , (2) 

n 

which expands the wavefunction as a linear combination 
of mutually orthogonal states 

=5^, (3) 

the states |$(Z)) which comprise the GCM are not mutu¬ 
ally orthogonal by design, namely the state overlap ma¬ 
trix 

N(Z';Z) = ($(Z / )|$(Z)), (4) 

is not diagonal. 

Of primary interest in this work are fermions (electrons 
or nucleons) for which each basis state |$(Z)) is then 
a Slater determinant, weighted by the complex valued 
scalar function /(Z) in Eq. Jl]). Thouless’ Theorem 
parametrizes these determinants by an exponential act¬ 
ing on a reference state |4> 0 ) 

|$(Z))=e 2 |$ 0 ), (5) 


where, schematically 

Z = Y,Z»%, (6) 

is a mono-excitation level generator with corresponding 
coefficients Z = (Z M ), is the product of two fermion 
creation and/or annihilation operators. The structure 
of Z and |$o) in turn define the space spanned by all 
possible GCM wavefunctions. If the generators Z are not 
sufficiently general, breaking one or more symmetries of 
the hamiltonian, then it is not possible to represent all 
symmetry preserving wavefunctions. 

Solutions for the GCM wavefunction are obtained in 
the time-independent case via the Schrodinger equation 

H I'I'gcm) = E |4>gcm) , (7) 

which can be solved via a matrix representation, resulting 
in the Hill-Wheeler-Griffin (HWG) equation 

J dZ H(Z'; Z) /(Z) = E j dZ N(Z'; Z) /(Z), ( 8 ) 

which is a generalized eigenvalue problem in continuous 
index form. Here 

H(Z , ;Z) = <$(Z , )|H|$(Z)) (9) 

N(Z';Z) = (<F(Z')|$(Z)). (10) 

A strength of the GCM representation is that evaluating 
the matrix elements Eq. ©, m, is inexpensive, even be¬ 
tween single Slater determinants corresponding to differ¬ 
ent basis sets or when also evaluating operators between 
such states, by using the Lowdin or Onishi formulae (see 
Ref. js] for the combined expression and original refer¬ 
ences). 

In practice, the generalized eigenvalue problem Eq. © 
can be solved by explicit discretization. In this case the 
matrix formulation is given by 
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H/=£N/ 


( 11 ) 
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for K sample points Z*- 1 ),..., Z^U and the, to be de¬ 
termined, weightings fk = /(Z^)). Using a naive dis¬ 
cretization, the computational demand scales exponen¬ 
tially with the dimension D : the number of generator co¬ 
ordinates Z. If one takes modestly TV « 10 points per di¬ 
mension then the number of samples is K = N D = 10 D , 
quickly becoming intractable. 

Applications in nuclear structure tend to use only a 
few ( D < 3 ) generator coordinates, selecting only 
those operators (such as the quadrupole operator) which 
correspond to physically motivated nuclear deformation 
modes in the expansion, rather than the full space of 
mono-excitation operators. This initial choice was due 
to the phenomenological origin of the GCM and appli¬ 
cations of interest. Subsequent investigations over the 
last sixty years have been hampered in no small part due 
to the poor scaling behaviour with the increase in the 
number of dimensions. This is also true even when the 
GCM formalism is only used as a means of symmetry 
projection for even a single broken symmetry Hartree- 
Fock (HF) or Hartree-Fock-Bogoliubov (HFB) state and 
the weightings /(Z) are fixed by symmetry. Other ap¬ 
proaches have parameterized the weighting function /(Z) 
by a combination of gaussian functions to restore some of 
the continuous flavour of the GCM ansatz, however the 
underlying numerical issues still remain. 

To circumvent these difficulties, we propose an alter¬ 
nate formulation of the GCM. This is possible by using a 
hole-particle representation of the generator coordinate 
operators which then allows for analytic symmetry pro¬ 
jection using tensor operators and the invariant mean 
(operator average) 11 E3- Symmetry projection then 
becomes a “back of an envelope calculation”. Explicitly, 
such mono-excitation generators Z of the hole-particle 
type are defined as 

Z = Y j Zai^i, ( 12 ) 

AI 

when generating the Hartree-Fock (HF) class of states 
!<F(Z)} and 

Z = 'y ' Zai a^a,/ 

AI 

+ \ ^ Z* AB a^ajg + \ ^ Z u a/aj , (13) 

AB IJ 

when generating the Hartree-Fock-Bogoliubov (HFB) 
class of states. Here I, J index the occupied spin orbitals 
(holes) and A , B the unoccupied spin orbitals (particles) 
with respect to the single-determinantal reference |<f>o)- 
The resulting GCM formulation has a direct connec¬ 
tion with projected Hartree-Fock (PHF) and projected 
Hartree-Fock-Bogoliubov (PHFB) wavefunctions, which 
then correspond to the single term approximation of the 
explicitly projected GCM. It is also related to the coupled 
cluster (CC) approach mm (see also Refs. EMI), a 
mapping from the general GCM prescription to CC can 


be constructed, even when no explicit symmetry break¬ 
ing is present. This mapping between GCM and CC 
leads to a distribution-like parameterization of the clus¬ 
ter amplitudes, and an alternate route to solve for the 
GCM wavefunction via the CC equations. The following 
sections will develop this new prescription. 

II. MAPPING FROM GCM TO CC 

In this section a mapping from the GCM representa¬ 
tion 

I'I'gcm) = y dZ/(Z) |$o), (14) 

to the coupled cluster representation 

|*cc> = e f |$ 0 >, (15) 

of the wavefunction will be elucidated. Here 

T = f 1 + T 2 + T 3 + T 4 + ... (16) 

are the conventional mono-, bi- and higher order cluster 
excitation operators (see e.g. Refs. [13, [lj])- Explicitly, 

exp(T) = 1 

+ tr 

+ T2 + ^Ti (17) 

+ T 3 + T 2 Ti + ^jTi 

+ T 4 + T 3 Tr + ^T| + ^T 2 T i + itj + ... , 

where terms corresponding to the same excitation order 
have been grouped. The mapping from GCM to CC can 
be expressed from two different perspectives on how the 
symmetry projection is accomplished: either when the 
weightings /(Z) are constrained via external means (such 
as via Lagrange multipliers), or when the weightings are 
unconstrained but an explicit projection has been per¬ 
formed, which the authors recommended as it is the most 
computationally efficient. These mappings from GCM to 
CC are also possible even when symmetry projection is 
not present or required. 

A. Constrained Weighting Formulation 

The GCM wavefunction Eq. ea is generally formu¬ 
lated in terms of broken symmetry generators Z as in 
Eqs. (TH21) or (fl3l) however the final wavefunction |II/gcm) 
is also implicitly assumed to have good quantum num¬ 
bers and therefore be symmetry adapted. Unless an ex¬ 
plicit projection has already been performed (as will be 
illustrated in the next section), then the weightings /(Z) 
must satisfy additional constraints to ensure symmetry 
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is preserved for the overall wavefunction. This can be 
accomplished by an implicit parameterization (if possi¬ 
ble) or by the use of additional constraints imposed via 
Lagrange multipliers. In either approach, the resulting 
GCM wavefunction is then given by 

I^gcm) = y dZ/ C (Z) |d> 0 ), (18) 

where f c ( Z) now denote the constrained weights which 
ensure symmetry is preserved. Possible methods to de¬ 
termine these weights include by explicit incorporation of 
the projection via a Baker-Campbell-Hausdorff (BCH) 
expansion (which also effectively modifies Z) or by re¬ 
quiring the symmetry violating contributions vanish at 
each order; similar to the Brillouin conditions for HFB. 
How these weights can be obtained will not however be 
considered further in this work. Instead a mapping to 
the coupled cluster form of the wavefunction will now be 
developed. 

The constrained weighting formulation of GCM 
Eq. m can be compared with the coupled cluster rep¬ 
resentation Eq. (ED, 

e f \<f> 0 )=JdZf c (Z)e i \$ 0 ), (19) 

for each excitation order Z fc , as compared to the corre¬ 
sponding fc-th row in the cluster expression in Eq. ED- 
At zeroth order, 

l = JdZf c (Z), (20) 

one obtains a normalization condition on the constrained 
weights f c (Z). One is tempted to treat / C (Z) as a distri¬ 
bution, however the scalar weights are themselves gener¬ 
ally complex-valued and even when real valued can have 
regions in the parameter space where the weight is neg¬ 
ative. 

Turning to the first order term in the expansion, 

f 1 = f dZf c (Z)Z, (21) 

this equality suggests the interpretation that under 
the GCM to CC mapping, the coupled cluster mono¬ 
excitation operator Ti corresponds to an average over 
the ‘pseudo-density’ / c (Z). Explicitly, in terms of the 
generator coordinate parameters, with Ti defined in a 
similar fashion to Eq. ©, 

Tji. = j dZi... dZ M ... dZn / C (Z) Z M 

= JdZpfMiZjZp, ( 22 ) 

where 

/ C M (ZJ = jdZx... dZ^-idZ^+i ...dZ D f c ( Z), (23) 


is a reduced (constrained) weight obtained by integrating 
over all generator coordinates Z except for Z From the 
normalization condition Eq. ED, then 

JdZ lt fM(Z li ) = 1, (24) 

for all choices of mono-excitation index /r. This again 
suggests the interpretation that / c (Z) can be considered 
a (complex valued) distribution, in which case /J^( Z M ) 
is then a one-parameter reduced weight. 

At quadratic (Z 2 ) order in the expansion, 

(f 2 + ±T^=±JdZf c (Z)Z 2 . (25) 

Rewriting 

T 2 = ^J dZf c (Z)7? ~ l -T{ 

= ±JdZf c (Z)(Z-T 1 ) 2 , (26) 

The bi-excitation operator T 2 obtained from the con¬ 
strained weighting GCM to CC mapping is the fluctua¬ 
tion (deviation) from the average. As before, one can for¬ 
mulate a two-parameter reduced weighting ZZ v ) 

for fj, 7 ^ v. 

Continuing in a similar fashion, one finds the higher 
order excitations 

T k = ^JdZf c (Z)(Z-f 1 ) k , (27) 

as the moments (or cumulants) in the space of mono¬ 
excitations about the mean given by Ti. This is the 
underlying structure of GCM, as cast in the constrained 
weighting formulation. 

B. Explicit Projection Formulation 

One difficulty with the approach taken in the previ¬ 
ous section is that the constrained weights f c ( Z) require 
either a specific parametrization or the use of external 
constraints, increasing the difficulty for a practical solu¬ 
tion. This constraint on the weightings can be removed if 
explicit projection is used. The weights /(Z) are then un¬ 
constrained and generally will not have the same values 
as f c { Z), although a mapping between these quantities 
is possible. 

Consider the explicit projection, 

P I^gcm) = P j dZ f(Z) |<E»o) 

= JdZ (p/(Z)e*|* 0 >). (28) 

If the GCM wavefunction is already symmetry adapted, 
such as by using constrained weights f c ( Z), then this is 
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identical to the original expression Eq. (fl4l) . If this is 
not the case, then evaluation of the projection even for a 
single state, 

P|$(Z))=Pe g |$ 0 ), (29) 

can be extremely computationally demanding or even in¬ 
tractable, depending on the symmetry group G of the 
projection operator P and which subgroup H C G the 
generators Z transform under (M3|. Spin angular mo¬ 
mentum G = SU(2) is the first example of such a ‘prob¬ 
lem group’, where exact state projection methods scale 
with the size of the state space, i.e. the number of possi¬ 
ble Slater determinantal configurations. 

If the reference |$o) is however already symmetry 
adapted and Z is of the hole-particle type (as we have 
explicitly chosen), then one can determine an analyti¬ 
cal form for the symmetry projection through the use of 
tensor operators and the invariant mean (operator aver¬ 
age) a lioj . Instead of the state projection or approxi¬ 
mate numerical projection, one employs operator projec¬ 
tion 

Usa = M G (e 2 ) , (30) 

in which case 

Pe 2 |$ 0 ) = U sa |$ 0 ). (31) 

Here A1 g(‘) is the invariant mean, which for discrete 
groups G corresponds to the average 

Xg(A) = ^ gAg” 1 . (32) 

11 gee 

For SU(2) spin projection, see Ref. 0- This equality 
Eq. (15T1) only holds under certain assumptions on the 
quantum numbers, for the case of SU(2) that the refer¬ 
ence is a closed-shell singlet 5 = 0, however this method 
can be extended to non-singlet quantum numbers sim¬ 
ply by inclusion of additional tensor operator coupling 
terms weighted with the appropriate Clebsch-Gordon co¬ 
efficient. The operator projection Eq. m then can be 
expanded for each order in Z, and significantly for each 
order the symmetry projection is algebraic and exact. 

The structure of the symmetry adapted correlation op¬ 
erator U sa can be formally linked to the coupled cluster 
expansion, so that 

U sa = exp(S), (33) 

where 

S = Si + S 2 + S 3 + S 4 + ... (34) 

are the symmetry adapted excitations, in correspondence 


with Eq. ©. Comparing Eqs. (1301) with (1331) . then 
S!=M g (Z) (35) 

S 2 =MgQz 2 ) -is? (36) 

S 3 =M G ^Z 3 j -S 2 S!-is? (37) 

S 4 = M g (iz 4 ) - S 3 Si - ^S 2 - is 2 S 2 - is 4 , 

(38) 

and so on. This approach was advocated in our previous 
work 0 as an exact method at any order to analytically 
determine projected broken symmetry wavefunctions, in¬ 
cluding for higher order operators approximations ( e.g. 
PUCCSD). Once the symmetry group G is known, ex¬ 
plicit forms for the invariant mean A4 G are relatively 
straightforward to evaluate for any fc-body operator, par¬ 
ticularly when Z is represented in terms of tensor opera¬ 
tors. 

As a simple example of this symmetry projection, con¬ 
sider the generator Z corresponding to real-valued BCS 
states, 

£(1V,S,M) 


N,S,M 

+1 


= £ zrA 

(39) 

N=—l 


given in terms of its tensor components Z^ N,S,M \ 
corresponding cluster operators S are given by 

The 

S 3 = z£'°’ 0) 

(40) 

A A(1,0,0) A( —1,0,0) 

b 2 — 

(41) 

S 3 =0 

(42) 

§4=-j (zr >0) ) 2 (zk- 1 '”- 0 ') 2 


_ 1 c2 

- 4 & 2» 

(43) 


which can be verified via Eqs. (Ha m- This can simi¬ 
larly be applied to other symmetry groups. Explicit ex¬ 
pressions for A4 g(-) operator projection of up to fourth 
order will be provided in a forthcoming paper fl 8 | for the 
symmetry group G = S x T x N, which corresponds to 
the projected real and complex-valued RHF, UHF, GHF, 
BCS and HFB classes of wavefunctions. 

For the remainder of this paper it is taken that the 
assumptions required for the invariant mean projection 
Eq. (1311) to be valid hold. For specificity, in case of the 
quantum chemistry literature this would require that the 
orbitals are real valued canonical spin-orbitals (space x 
spin) and that the reference |$ 0 ) is a closed shell sin¬ 
glet state with good particle number. For simplicity, the 
generator coordinates Z and weightings /(Z) will also be 
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assumed real valued. These conditions are intended to 
simplify the derivation of this alternate GCM prescrip¬ 
tion, avoiding the potential quagmire with spin coupling 
and other issues for which a unique definition of the hole- 
particle representation is not as straightforward. Some of 
these conditions can be relaxed and the GCM structural 
form generalized further, such as by allowing for different 
reference states or by choosing |4>(Z)) to include both sin¬ 
gle and double excitations for the generator coordinates. 
The authors however wish to first grab the proverbial 
low-lying fruit first. 

Using the operator average, then 

P|*gcm> = P fdZ /(Z) M G (e 2 ) |$„>, (44) 

or using Eqs. (I5T1) . (1351) . 

P|*gcm> = J dZ /(Z) e § ^ |*o> ■ (45) 

In this representation, the GCM is a (continuous) linear 
combination of cluster operators, each corresponding to 
a projected broken symmetry wavefunction. 

As before, comparing order by order with the coupled 
cluster expansion 

l=JdZf(Z), (46) 

but now for the unconstrained weightings. The mono¬ 

excitation operator is given by 

T i = J dZf(Z)M G (Z) 

= y'dZ/(Z)S 1 (Z), (47) 

and so the unconstrained mono-excitation Ti is an aver¬ 
age over projected mono-excitation generator Z. 

Turning to the bi-excitations, 

= dZ /(Z)M G (Z 2 ), (48) 

or 

T2 = 1/ dZ/(Z)X G (Z 2 )-it?. (49) 

The bi-excitation can be interpreted, just as before, as 
the deviation from the average amplitude 

T 2 = i J dZ /( Z) M g ((Z - Ti) 2 ), (50) 

where however the symmetry projection via the invariant 
mean A4 G (-) has been incorporated with unconstrained 
weightings /(Z). 


Alternately, the bi-excitations can be re-cast using the 
cluster expansion for each term, 

T2 = | dZf(Z) Qs 2 (Z) + S 2 (Z)) - if 2 
= j dZf(Z) S 2 (Z) 

+ ±JdZf(Z)(s 1 (Z)-T i y, (51) 

in which case the T 2 comes both from the sum of a 
weighted S 2 and the deviation from the operator average 
of Ti. Here the S 2 contribution comes from explicit sym¬ 
metry breaking, whereas the second term comes from the 
spread of the mono-excitation coefficients or alternately 
the cross-correlation. The S 2 term vanishes entirely if no 
symmetry breaking is present and can be interpreted as 
the canonical order parameter. One finds that T 2 only 
vanishes if both S 2 = 0 and 

/(Z) = S(Z — Z 0 ) , (52) 

that is that the GCM corresponds to a single symme¬ 
try restricted Slater determinant, i.e. restricted Hartree- 
Fock (RHF) wavefunction. 

At the T 2 bi-excitation level, GCM is a (continuous) 
linear combination of projected broken symmetry wave- 
functions as in Eq. (1451) or instead via Eq. ED is a 
weighted distribution of pairing amplitudes. The latter 
scheme can be related to the plethora of approximate 
doubles schemes, such as perfect pairing, RVB, GVB, 
canonical BCS, AGP, and so on. This representation 
thus provides a natural framework to discuss the spar¬ 
sity structure of the cluster amplitudes, which can be 
dictated by symmetry or other considerations. 

For the third order expression, 

T 3 + t 2 T! + It 3 = I J dZ /(Z)A4 G (Z 3 ). ( 53 ) 

Re-arranging, 

T3 = |/ dZ /(Z) A4 g ((Z - TO 3 ), (54) 

which can be linked directly to the corresponding symme¬ 
try constrained weighting / C (Z) expression. Alternately, 
in terms of the symmetry adapted cluster operators, 

T 3 = JdZf(Z)S 3 (Z) 

+1 dZf(Z) (S 2 (Z)-T 2 )(S 1 (Z)-T 1 ) 

+ I|dZ/(Z)(S 1 (Z)-T 1 ) 3 . (55) 

The expansion suggests that the higher order cluster am¬ 
plitudes, when interpreted through the GCM to CC map¬ 
ping, are a combination of ‘pure’ symmetry breaking con¬ 
tributions together with distributional (cross-correlation) 
contributions. 
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Fourth order proceeds similarly, 

T 4 = i I dZ /(Z) M g ((Z - Ti) 1 2 3 4 5 ), (56) 

and for the cluster formulation 

T 4 = /dZ/(Z)S 4 (Z) 

+ j dZ f(Z) (S 3 (Z)-T 3 )(Si(Z)-f 1 ) 

+ 4/ rfZ/(Z)(S 2 (Z)-T 2 ) 2 

+ 4/ rfZ/(Z)(S 2 (Z)-T 2 )(S 1 (Z)-T 1 ) 2 

+ l|dZ/(Z)(S 1 (Z)-T 1 ) 4 (57) 

This is the structure of the GCM under the explicitly pro¬ 
jected CC mapping, which provides an alternate route to 
a highly efficient and direct means of symmetry projec¬ 
tion. 


III. FIRST CHECKPOINT 

This work has provided a clear prescription for how 
to utilize the structural form of GCM while avoiding the 
historical computational bottleneck. This has been ac¬ 
complished through the use of analytical projection via 
the invariant mean and tensor operator representations, 
which can be effected efficiently when the generator co¬ 
ordinate is cast within the hole-particle representation. 

At the lowest level of approximation, which takes a 
single Slater determinant for the GCM, this approach 
yields an efficient method to obtain the wavefunction for 
individual PHF and PHFB states as cast within the cou¬ 
pled cluster representation. The imposed sparsity on the 
cluster amplitudes (for a single or linear combination of 
Slater determinants) can be formulated in this frame¬ 
work, and related to the canonical (finite size) order pa¬ 
rameters of the system. For higher order amplitudes, 
the symmetry projection can also be used to obtain the 
externally corrected correlation corrections for the T 3 


and T 4 amplitudes. These can then be determined self- 
consistently through the higher order cluster equations, 
as significantly fewer free parameters are required as than 
a full T 3 or T 4 calculation. 

The GCM expansion provides a natural approach to 
go beyond projected HF and projected HFB wavefunc- 
tions towards the exact solution, and in this work has 
been mapped to the coupled cluster amplitudes. Un¬ 
der this mapping, the resulting GCM weights /(Z) can 
be interpreted as a pseudo-distribution. Provided the 
weights /(Z) between uncoupled subsystems are uncorre¬ 
lated, then the GCM wavefunction is size-extensive, just 
as for the coupled cluster approach. 

An explicit solution procedure has not been presented 
in this work, however one can use the coupled cluster 
equations directly as the authors have advocated previ¬ 
ously for projected broken symmetry wavefunctions 11- 
Even at the single term approximation of GCM, which 
corresponds to the PBF or PBFB wavefunctions, this 
leads to a reduction of the cluster equations to a rela¬ 
tively few expressions for the free parameters. This ap¬ 
proach can also be used for GCM even when symmetry 
breaking is not present, such as when choosing D physi¬ 
cally motivated generator coordinates. The coupled clus¬ 
ter equations then give a D 6 7 8 scaling, once intermediates 
have been calculated. 

The mapping from GCM to CC provides one notewor¬ 
thy point for the coupled cluster approach, illustrating 
the severity of conventional coupled cluster truncation 
schemes. In the case of CCSD, singles and doubles am¬ 
plitudes Ti and T 2 are evaluated but all higher order 
amplitudes Tfc for k > 2 are omitted (zeroed). This 
places a strong constraint on the pseudo-distribution 
/(Z) which results in unphysical artifacts in the resulting 
self-consistent CCSD equations. One avenue for further 
investigation is how the pseudo-distributional structures 
for T 3 and T 4 can be exploited to address this deficiency. 
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